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ABSTRACT 

Context. The evolution of a circumstellar disk from its gas-rich protoplanetary stage to its gas-poor debris stage is not understood well. 
It is apparent that disk clearing progresses from the inside-out on a short time scale and models of photoevaporation are frequently 
used to explain this. However, the photoevaporation rates predicted by recent models differ by up to two orders of magnitude, resulting 
in uncertain time scales for the final stages of disk clearing. 

Aims. Photoevaporation theories predict that the final stages of disk-clearing progress in objects that have ceased accretion but still 
posses considerable material at radii far from the star. Weak-line T Tauri stars (WTTS) with infrared emission in excess of what 
is expected from the stellar photosphere are likely in this configuration. We aim to provide observational constraints on theories of 
disk-clearing by measuring the dust masses and CO content of a sample of young (1.8-26.3 Myr) WTTS. 

Methods. We used ALMA Band 6 to obtain continuum and 12 CO(2-l) line fluxes for a sample of 24 WTTS stars with known infrared 
excess. For these WTTS, we inferred the dust mass from the continuum observations and derived disk luminosities and ages to allow 
comparison with previously detected WTTS. 

Results. We detect continuum emission in only four of 24 WTTS, and no 12 CO(2-l) emission in any of them. For those WTTS where 
no continuum was detected, their ages and derived upper limits suggest they are debris disks, which makes them some of the youngest 
debris disks known. Of those where continuum was detected, three are possible photoevaporating disks, although the lack of CO 
detection suggests a severely reduced gas-to-dust ratio. 

Conclusions. The low fraction of continuum detections implies that, once accretion onto the star stops, the clearing of the majority of 
dust progresses very rapidly. Most WTTS with infrared excess are likely not in transition but are instead young debris disks, whose 
dust is either primordial and has survived disk-clearing, or is of second-generation origin. In the latter case, the presence of giant 
planets within these WTTS might be the cause. 

Key words. Protoplanetary disks - Planets and satellites: formation - Planet-disk interactions - Radio continuum: planetary systems 
- Radio lines: planetary systems - Infrared: planetary systems 


1. Introduction 


Circumstellar disks typically fall into two categories - the mas¬ 
sive gas-rich protoplanetary disks and the gas-poor debris disks. 
Protoplanetary disks are a natural consequence of the star for¬ 
mation process, with the majority of their evolution being dom¬ 
inated by viscous accretion onto the star (e.g., |Lynden-Bell &| 
Pringle 1974). They can easily be identified by their near and 
mid-infrared (mid-IR) excesses, caused by their vast optically 
thick dust disks (e.g. Strom et al. |1989) >, although the majority 
of their mass in fact resides in gas with typical gas-to-dust mass 
ratios of ~100 (e.g. |Pietu et al.|2005[|Panic et al.|2008] l. The de¬ 
bris disks, meanwhile, generally contain little or no detectable 
gas, with dust often confined to a narrow ring (see|Wyatt|2008 


for a review). It has long been suggested that debris disks could 


be a later stage of evolution from the protoplanetary disks, but 
the nature of the main physical processes that drive this evolu¬ 
tion is ill-understood and remains one of the biggest questions in 
the field. To understand this phenomenon, the nature of a class 
of disks known as transition disks, needs to be investigated. 

While a unique and universally accepted definition of what 
constitutes a transition disk (TD) does not exist, the most general 
definition is that of a T-Tauri disk with reduced excess emission 
at near to far IR wavelengths relative to typical T-Tauri disks. 
Most TDs show significantly reduced fluxes at short wavelengths 
(< 10 pm), while still showing average emission at longer wave¬ 
lengths. This observation is indicative of a dust cavity in the in¬ 
nermost region of the disk, and indeed sub-mm images of TDs 
have proved this to be the case ( |Pietu et al.|2006||Hug hes et al. 
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|2007||Brown et ak]2 009; Ca novas et al.|2015) . This kind of ge¬ 
ometry requires a mechanism that can clear the dust from the 
inside out, which is not consistent with the previously domi¬ 
nant mechanism of viscous accretion (Hart mann et al.|1998[|Ar-| 
mitage et al. 1999|l. In addition, the lifetime of this transition 


stage must be short, i.e. <1 Myr because the detection rate of 
TDs is comparatively low, which again argues against the slow 
process of viscous accretion ( Wolk & Walter) 1996) ; Duvert et al. 


2000; Andrews & Williams 2005 2007| l. Many mechanisms 
have been proposed to explain the phenomena of TDs, including 
grain growth, giant planet formation, binarity and photoevapora¬ 
tion. Although all these mechanisms probably contribute to disk 
evolution, photoevaporation is the most plausible mechanism to 
explain the final rapid removal of material at large radii. 


Initial theories of photoevaporation (e.g. Hollenbach et al. 
|1994| |Clarke et al.||2001] l described the photoionisation of hy¬ 
drogen in the disk surface by extreme ultraviolet (EUV) pho¬ 
tons. This photoionisation forms a pressure gradient that is able 
to drive mass loss in a photoevaporative wind beyond a critical 
radius. Evidence to support this theory has been found by com- 


paring models to observations of both [Nell] ( 

Alexander||2008; 

Pascucci & Sterzik 2009;|Ercolano & Owen 20 

0; Pascucci et al. 

2011 

Sacco et al. [2012) and [OI] lines (Font et al. 2004 Gorti 

et al. 

2011 

Rigliaco et al. 2013)). Although this process is likely 


to occur throughout the disk’s lifetime, it only becomes signifi¬ 
cant when the accretion rate becomes comparable to the photoe¬ 
vaporation rate. When this happens, photoevaporation can open 
a gap, forming an inner and an outer disk. The inner disk, now 
cut-off from re-supply, drains on a viscous time scale, thus cre¬ 
ating the observed transition disk geometry (see Alexander et al. 
|2014| for a review). In the case of a transition disk however, the 
inner opacity hole allows radiation from the star to then reach 
the rim of the outer disk unimpeded, allowing it to complete the 
clearing of the disk on a short timescale of approximately 10 5 


years (Alexander et al. |2006[ Alexander & Armitag e|2007) l. 

Although a very promising theory, EUV photoevaporation 
suffers from an uncertainty in the EUV flux incident on the disk. 
Stellar EUV flux is difficult to measure as interstellar absorp¬ 
tion prohibits direct observation of the ionising photons. Further¬ 
more, these ionising photons can be blocked by optically thick 
accretion columns, jets or winds, resulting in final estimates for 
the flux that can reach the disk varying by orders of magnitude 
( |Herc zeg 2007 |Pascucci et al.|20l2) . Some recent advances in 
this area have had some success using free-free emission from 
the disk to place limits on the EUV flux, but the results suggest 
EUV wavelengths are, in fact, not sufficient to explain the [Nell] 
emission seen in some systems (Pascucci et al .|2014 1. 

More recently, the effect of both X-ray (Owen et al. pOlO 


2011 


2008 


2012} and far ultraviolet (FUV) ( |Gorti & Hollenbach | 


2009 Gorti et aI7||2009 1 wavelengths have been included 


into photoevaporation models and a lot of debate surrounds the 
question of which wavelength is most dominant. X-ray and FUV 
both predict the same evolutionary behaviour as EUV, but gener¬ 
ally predict mass loss rates which are orders of magnitude higher 
than the original EUV models, with M PE ~ 10 1,1 M Q yr for 
EUV models, and M PE ~ 10 8 M G yr '' for FUV/X-ray (Gorti 


et al. |2009[ ). Although stellar FUV and X-ray emissions are eas¬ 
ier to measure than EUV, these models also suffer from uncer¬ 
tainties regarding disk chemistry and dust properties. FUV mod¬ 
els, for example, describe heating of the gas being dominated 
by polycyclic aromatic hydrocarbons (PAHs), but the abundance 


and depletion of PAHs are difficult to determine (Geers et al. 
12009) ). 
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To answer some of the questions surrounding photoevapo¬ 
ration models, a study of transition disks in their last stages of 
gas clearing is required. A natural sample in which to find such 
disks is in weak-lined T-Tauri stars (WTTS). Unlike the classical 
T-Tauri stars (CTTS), WTTS have a narrow Her width, which is 
a strong indicator that the star is no longer accreting (or is ac¬ 
creting only at a relatively low level), and therefore lacks gas at 
radii close to the star (see also Pascucci et al. 2006; Ingleby et al. 


2009). The majority of WTTS also have no IR excess, suggest¬ 


ing they have already cleared all their circumstellar gas and dust 
and are not in transition. Nevertheless, a relatively small percent¬ 
age of WTTS (~ 20%) do display an IR excess, and this excess 
suggests a significant amount of dust (Cieza et al. 2007; Wah- 
haj et al. 2010). Even so, a measurement of their dust mass has 
not yet been achieved and their nature as transition disks has not 
been confirmed. 

For example, in a dedicated survey of transition disks se¬ 
lected from Spitzer ( |Cieza et ah 2010 2012| Romero et al.|2012 1 
only one such WTTS system was detected at sub-mm wave¬ 
lengths (FW Tau), but observations of this particular system with 
ALMA recently revealed that the sub-mm emission originates 
from around an accreting third object. FW Tau should therefore 


not be considered as a disk-possessing WTTS (Kraus et al. 2015 


Caceres et al. 2015). Furthermore, no previous studies have been 
able to measure the gas content of these disks, and it is still to 
be determined whether WTTS with an IR excess (referred to 
as IR-WTTS in the remainder of this paper) are still in transi¬ 
tion and photoevaporating their gas, or whether they are more 
akin to young debris disks. If they should come under the former 
class, these systems will indeed be in the final stages of gas clear¬ 
ing, and observational constraints on their mass and gas-to-dust 
mass ratio will be invaluable to theories of photoevaporation. If 
they should come under the latter class, then these disks will be 
among some of the youngest debris disks known, and their ages 
can provide a constraint on the initial conditions of the debris- 
disk phenomenon. Furthermore, if this result were to hold for a 
large sample of IR-WTTS, then this would suggest that by the 
time accretion ceases, the gas has either already been heavily 
depleted or the remaining gas has been photoevaporated rapidly. 
In this sense, such a result would have a strong influence on both 
photoevaporation theory and disk evolution in general. 

Here, we present one of the first detailed studies of WTTS 
with IR excesses using ALMA. We use Band 6, observing both 
the continuum and the 12 CO(2-l) transition with the aim of iden¬ 
tifying the evolutionary state of these WTTS disks and imposing 
limits on photoevaporation theory. 


2. Observational procedure 

Characterisation of the disks around IR-WTTS can be achieved 
by measuring their dust and gas mass, stellar age, fractional disk 
luminosity and multi-wavelength photometry. It is these values 
that we therefore attempt to determine. 


2.1. Target selection 

Our sample consists of 24 pre-main-sequence stars (based on 
weak Her emission and Li I absorption) in nearby (< 200 pc) 
molecular clouds (see Table [lj. The objects have been classed 
as a WTTS based on the velocity widths of their Her lines. The 
width of the Her line provides a reliable, distant-independent in¬ 
dication of a star’s accretion, with accretion producing broad, 
asymmetrical Her emission. Non-accreting objects will still pro¬ 
duce Her emission of chromospheric origin, but this tends to 
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Table 1 . Sample parameters 


No. 

2MASS ID 

Cloud 

Distance 

(parsecs) 

Spectral type 

H„ width 
(kms' 1 ) 

Ref. 

Binary sep. 
(arsecs) 

Binary ref. 

1 

04182147+1658470 

Taurus 

135+20 

K5 

Absorp. 

1,2,3 



2 

04192625+2826142 

Taurus 

135+20 

K7 

180 

4,5 

10.5 

1 

3 

04242321+2650084 

Taurus 

135+20 

M2 

200 

5 



4 

04314503+2859081 

Taurus 

135+20 

F5 

Absorp. 

5 



5 

04325323+1735337 

Taurus 

135+20 

M2 

138 

4,6 



6 

04330422+2921499 

Taurus 

135+20 

B9 

Absorp. 

5 



7 

04364912+2412588 

Taurus 

135+20 

F2 

Absorp. 

2,5 



8 

04403979+2519061 

Taurus 

135+20 

M5 

130 

3,5 



9 

04420548+2522562 

Taurus 

135+20 

K7 

94 

3,6 

0.3 

11 

10 

08413703-7903304 

7/Cham aleonis 

97+3 

M3 

Absorp. 

4,7 



11 

08422372-7904030 

7/Cham aleonis 

97+3 

M3 

Absorp. 

4,7 



12 

11073519-7734493 

Chamaleon I 

160+15 

M4 

89 

1,8 



13 

11124268-7722230 

Chamaleon I 

160+15 

G8 

Absorp. 

1,8,9 

0.247 

12 

14 

16002612-4153553 

Lupus IV 

150+20 

M5.25 

162 

10 

2.8 

13 

15 

16010896-3320141 

Lupus I 

150+20 

G8 

Absorp. 

6 



16 

16031181-3239202 

Lupus I 

150+20 

K7 

132 

4,6 



17 

16085553-3902339 

Lup III 

200+20 

M6 

189 

10 

2.8 

13 

18 

16124119-1924182 

Ophiuchus 

119+6 

K8 

131 

6 

1.0 

14 

19 

16220961-1953005 

Ophiuchus 

119+6 

M3.7 

132 

10 

1.8,3 

10 

20 

16223757-2345508 

Ophiuchus 

119+6 

M2.5 

128 

4,6 



21 

16251469-2456069 

Ophiuchus 

119+6 

M0 

206 

4,6 



22 

16275209-2440503 

Ophiuchus 

119+6 

K7 

151 

4 

0.480 

15 

23 

19002906-3656036 

Corona Australis 

129+11 

M4 

93 

10 

0.132 

16 

24 

19012901-3701484 

Corona Australis 

129+11 

M3.75 

83 

10 

0.5 

10 


Notes. The target’s photometry, spectral types and H tt widths can be found in the following papers:- l)|Nguyen et al.| |2012|); 2) |Howard et al. 
j~3); 3)|Luhman et al.|p0T0l>; 4 ) [CiezaetaT|(|20T3[; 5)|Cieza et al.|pOT2> ; 6) |Wahhaj et aflpUTO) ; 7) |Sicilia-Aguilar et al.|j2009>; 8) '~ ' 


Luhman 


|et al.jl 20081; 9) Matraet al. (2012); 10) Romero et al.| 201 2}. The binary systems were identified in the following papers:- 1 l)|Leinert et a .|(|1993|t; 
12) ILafreniere et al.| ( |2008l >; 13) |Merin et al.| ( |2008| ); 14) |Prato| ( |2007> ; 15) |Ratzka et aT71 ( [2005l >; 16) |Kohler etlflpMj . 


be comparatively narrow and symmetric. The empirical divid¬ 
ing line between accreting and non-accreting has been a mat¬ 
ter of some debate, with some claiming thatjiccreting systems 
have a 10% peak width of >270 kms -1 (White & B asri[2 003l>, 
and others claiming that that the dividing line is >200 kms" , 
which varies with spectral type ( |Martm|jl9981 >. The systems in 
this study were selected because they all either have 10% peak 
widths less than 200 kms -1 , or instead display Her absorption 
lines, allowing them to be identified as very likely non-accretors. 
They all lack considerable IR excess at ~10 gm or shorter wave¬ 
lengths but show weak yet robust (>5-10 cr) excesses in the mid 
and/or far-IR from Spitzer, WISE, and/or Herschel. These targets 
have been previously labelled as either photoevaporating transi¬ 
tion disks, or debris disks candidates, because of their lack of ac¬ 
cretion and because of their estimated fractional disk luminosi¬ 
ties having value s less than the typical value for protoplanetary 


disks of ~ 10 (Wahhaj et al. 


2010 


Cieza et al. 


2010 


2012 


Rom ero et al.|2012[ ). Those with fractional disk luminosities in 
the region of 10 1 > Ld / L, > 10 2 were classed as photoe¬ 
vaporating transition disks, while those with a lower value were 
classed as debris disk candidates. Binaries with wide projected 
separations were included, but no spectroscopic or confirmed 


close binaries were, because circumbinary disks are believed to 
undergo a different evolution and therefore need to be studied 
separately (Kraus et al .|2012| l. System 23 is a possible close bi¬ 
nary, but only single epoch data is available for this system and 
its binary nature has yet to be confirmed via proper motion. 

2.2. ALMA observations 

Observations of the above systems were performed in Band 6, 
with Systems 14, 17, 19, 23 and 24 being observed in Cycle 0 
(2012), and the remaining in Cycle 1 (2013). Cycle 1 observa¬ 
tions were split into three, based on their host cloud, with the Lu¬ 
pus and Ophiucus systems in one group, Taurus in another, and 
Chamaleonis in the final group. These groups were observed at 
different times and with slightly different configurations. Table[2] 
gives a summary of these observational set-ups. 

We chose to observe the 12 CO(2-l) line since it is highly 
sensitive to the presence of circumstellar gas out to large radii, 
where the bulk of the gas should be located. We obtained one 
epoch of observation for all systems, with the correlator config¬ 
ured to obtain one baseband centred on 230.52 GHz which was 
aimed at detecting the 12 CO(2-l) spectral line, and three con- 
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Table 2. Observation log 


Date 

Cloud 

Antennas 

ToT 

(min) 

Antennas 

flagged 

BC 

(QSO) 

GC 

PC 

(QSO) 

PWV 

Min/Max 
baseline (m) 

12/06/17 

Lup/Op/CrA 

20 

4.54 


J1733-1304 

Neptune 

J1924-2914 

1.41 

21.2/402 

13/11/01 

Lup/Op 

29 

2.07 

DV19,08,06 

J1924-2914 

Neptune 

J1625-2527 

1.60 

17.3/992 

13/12/04 

Taurus 

27 

3.60 

DV08 

J0423-0120 

J0510+180 

J0510+180 

5.11 

15.8/463 

13/12/18 

Cha 

23 

4.16 

DV 19,08 

J1107-4449 

Ceres 

J0635-7516 

2.31 

15.1/992 


Notes. ToT - Time on target, BC - Bandpass calibrator, GC - Gain calibrator, PC - Phase calibrator, PWV - Precipitable water vapour 


tinuum basebands centred at 228.52, 214.52, and 212.52 GHz. 
However, a fault with the local oscillator during the Cycle 0 
observations meant that only the 230.52 and 228.52 GHz base¬ 
bands could be observed. The total bandwidth for the observa¬ 
tions was 3.75 GHz for the Cycle 0 and 7.5 GHz for Cycle 1 
observations, with a unique spectral resolution of 976.56 kHz in 
3840 channels for each 1.875 GHz baseband. In all cases, the re¬ 
quested rms was set at 0.16 mJy for the continuum, and 30 mJy 
for each individual velocity channel (and thus for the 12 CO(2- 
1) line). Standard calibration steps were applied to the data, and 
we obtained the final images by deconvolving the set of visibili¬ 
ties with the clean task implemented in CASA (McMullin et al. 
2007[>, using natural weighting. A point-source, fitted to the mea¬ 


sured visibilities, was used for those systems with a continuum 
detection to gain accurate estimation of this flux value. 


3. Results 

3.1. Disk dust masses 


Only four of our 24 systems have detectable continuum above 
the 3<x level, with an additional tenuous detection in system 8 . 
This system displays weak emission, but which is just below the 
3 sigma level and therefore not significant enough to claim a 
detection. 

Andrews & Williams (2005)> show that, owing to continuum 


emission being optically thin at mm wavelengths, the dust mass 
can be estimated by a simple equation of the form Mj ust — C v x 
F v where C v is a constant for a given frequency F v . We adopt the 
constant derived for 1.3mm by Cieza et al.] ( j2008| ), and use the 
equation 


M dust = 0.566 x 


F v ( 1300) 

f d fl 

mJy 

[UOpcJ 




( 1 ) 


to estimate the dust mass for our targets, or the dust mass upper 
limit in the case of a non-detection. This approach has signifi¬ 
cant uncertainties but is quite standard in the field and therefore 
allows for meaningful comparisons with previous results. The 
results of this are displayed in Table [3] In the case of a non¬ 
detection, we quote the 3 sigma value as our upper limit, with 
this uncertainty being dominated by the uncertainty in the dis¬ 
tance to the system. 


Table 3. Calculated parameters 


No. 

1.3mm Flux 

Dust mass 

Stellar mass 

Age 

Ld ! L* 


(mJy) 

(M e ) 

(Mo) 

(Myr) 

(xlO -3 ) 

1 

<0.436 

<0.30 



<3.6 

2 

0.533 

0.28+0.09 

0.67 

o c+2.7 
^ -1.1 

2.3 

3 

<0.426 

<0.30 

0.35 

6 7 + 2 - 7 

<1.5 

4 

<0.435 

<0.30 



<7.1 

5 

<0.438 

<0.30 

0.36 

3 1 +1 

<3.5 

6 

<0.431 

<0.30 



<2.9 

7 

<0.435 

<0.30 



6.5 

8 

<0.431 

<0.30 

0.17 

4 l +lA 

^ -0.7 

<11 

9 

<0.423 

<0.29 

0.65 

1 0+1.2 
-0.5 

<4.6 

10 

<0.476 

<0.14 

0.27 

9 0 +0 - 7 

7 -0.6 

<1.4 

11 

<0.472 

<0.14 

0.31 

o 4+0.5 
^ -0.4 

< 0.8 

12 

<0.541 

<0.48 

0.195 

15 . 0/3 

<2.3 

13 

<0.494 

<0.44 

1.4 

12 4 +2 
- 1.6 

<2.3* 

14 

0.696 

0.45+0.13 

0.17 

6.9/™ 

80 

15 

<0.442 

<0.37 

1.1 

26.3/™ 

<0.26 

16 

<0.453 

<0.38 

0.70 

7 7 + 2 - 3 
'■'-1.6 

< 2.8 

17 

1.813 

2.09+0.44 

0.10 

2 6 +1 - 2 

140 

18 

<0.459 

< 0.21 

0.64 

6 6 +2 - 4 
u.u 1 5 

<3.4 

19 

<0.483 

< 0.22 

0.34 

9 9+0.3 
^ -0.2 

<5.95 

20 

<0.461 

< 0.21 

0.33 

q q+0.8 
y - y -0.6 

<1.96 

21 

<0.453 

< 0.20 

0.56 

3 0 +0 - 4 

J - 0.6 

< 2.0 

22 

<0.442 

< 0.20 

0.77 

4 2 +L2 
- 0.8 

<1.5 

23 

0.569 

0.27+0.05 

0.24 

4 7+1.0 

^-0.7 

28 

24 

<0.465 

<0.26 

0.21 

10 . 2 /™ 

<5.7 


Notes. Four systems have no calculated stellar mass or age as their lu¬ 
minosity does not coincide with any PMS tracks and are therefore prob¬ 
ably background main-sequence stars. The final column of L 0 /L, gives 
the value of fractional disk luminosity, calculated by integration of the 
system spectral energy distribution. ‘The value of L c /L, < 1(T 3 quoted 
for system 13 assumes the far-IR emission is due to contamination. 


3.2. i2 CO(2-1) non-detections 

We did not detect 12 CO(2-l) emission for any of the observed 
systems suggesting that very little gas remains in the disks. De¬ 
termining upper limits on the gas mass that these non-detections 
imply is more difficult than for the dust, however, because the 
CO emission is generally optically thick and the conversion to 


H 2 (the component that makes up the majority of the gas mass) 
is uncertain. This uncertainty is largely due to photo-dissociation 
and freeze-out, which depend on the density and temperature 
structure of the disk. To help constrain disk gas masses from CO 
observations, Williams & Best ( 2014| ) created a grid of models 
with different density structures that include a simple prescrip¬ 
tion for the CO chemistry. The best constraints on gas masses 
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Log wavelength (/zm) 


Fig. 1. SEDs of the 4 targets in which 1.3mm flux was detected. The filled circles are archival photometry points, with the upper limits as arrows. 
The blue diamonds show the flux before being corrected for extinction. The solid green lines represent the stellar photosphere normalised to the 
extinction-corrected J band. The dotted lines correspond to the median mid-IR SED of K5-M2 CTTSs calculated by jFurlan et al.| ( f2006| >, with the 
dashed lines as the upper and lower quartiles. The red line corresponds to IRS spectra, where available. 


come from a combination of the moderately optically thin CO 
isotopologues but, because of the sensitivity of these observa¬ 
tions, we can obtain useful limits just by using the 12 CO(2-l) 
line. We compare our results to this grid for a stellar mass of 0.5 
M g , as this mass is suitable for the stars in our sample (see Table 
[3J, but we leave all other parameters free (disk mass, disk ra¬ 
dius, radial power-law index, surface-density power-law index, 
temperature profile, and inclination). For the most massive dust 
disk in our sample, i.e. System 17, an ISM gas-to-dust ratio of 
100 would imply a gas mass of ~ 1 Mj, which would be read¬ 
ily detectable in l2 CO(2-l) at our 3cr detection limit of 90 mJy 
kms~'. We can therefore conclude that the gas-to-dust ratio is 
substantially lower than 100 in this disk. A similar argument ap¬ 
plies to system 14. For the other disks, despite the enhanced ef¬ 
fects of photo-dissociation resulting from their lower maximum 
gas masses of M gas ~ 0.1 Mj, the CO emission would still be 
above our detection threshold for 90% of the models. This sug¬ 
gests that the gas-to-dust ratio in these disks has also evolved to 
well below the ISM value and that the amount of gas in these 
disks is extremely small. 

3.3. Stellar ages 

To be able to make a comparison with previously studied sys¬ 
tems, and to impose limits on disk-clearing timescales, the stel¬ 
lar ages are required. These were obtained by calculating stellar 


luminosities and temperatures, and comparing these to pre-main 
sequence (PMS) isochrones. Stellar temperatures were estimated 
by their spectral type, using the scale provided by Ken yon &| 
Hartmann ( fl995j >. The stellar luminosity was calculated by first 
applying a de-reddening correction to each star, as appropriate 
for its spectral type, again in accordance with the values pro¬ 
vided by Kenyon & Hartmann) ( |1995[ ). The J band magnitude 
was used as a reference because this band is less affected by ex¬ 
tinction than at shorter wavelengths, whilst having little chance 
of being affected by the flux from any circumstellar material. The 
distances in Table[l]were then used to calculate the stellar radius 
required to recreate the measured emission, and these radii used 
to calculate the stellar luminosity. The values of temperature and 


luminosity were then compared to the PMS isochrones of Siess 
|et al.[p000| to obtain the age. The results of this process are 
displayed in Table [3] along with the corresponding stellar mass 
suggested by the PMS isochrones. Four systems (1,4, 6, and 7) 
appeared too under-luminous for their temperature to coincide 
with any PMS model, and therefore do not have a corresponding 
value for the age column in Table [T] 


Some caution must be taken when interpreting these ages, 
since the evolutionary tracks for these PMS systems lie close to¬ 
gether on the Hertzsprung-Russel diagram. Moreover, the dis¬ 
tances to the individual objects have an uncertainty of up to 
—15%, which can introduce a ~30% uncertainty into the intrinsic 
luminosity of the object. These effects can conspire to result in a 
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high degree of uncertainty in the ages of individual stars. Fortu¬ 
nately, however, the ages of WTTS are considerably easier to de¬ 
termine than those of CTTS. This is because CTTS are affected 
by veiling and possess highly heterogeneous photospheres due 
to their accretion, which makes their intrinsic luminosities and 
temperature much more challenging to determine (Cieza et al. 
|2005) >. In the case of our WTTS, the largest contributor to the 
uncertainty is probably the distance of the object. 


3.4. Fractional disk luminosities 


A frequently used criterion for distinguishing between proto¬ 
planetary and debris disks is the fractional disk luminosity, 
Ld / L ,, which measures the fraction of the stellar radiation that 
is intercepted and re-emitted by the disk. Typical values for pro¬ 
toplanetary disks are in the region Lo/L* ~ 0.1 (e.g. Cieza et al. 
2010), whereas the value for deb ris disks is much l ower, with 
typical values Ld/L* < 10 3 (e.g. Decin et al7||2003 I. To deter¬ 
mine this quantity for our sample, we attribute all emission from 
the star-subtracted SED to the disk (Ld) and fit this emission to 
either one or two black bodies, as required. The black body emis¬ 
sion at wavelengths longer than To was modified by a factor of 
Uo/Af, since, at these wavelengths, the emission from an opti¬ 
cally thin debris disk is observed to deviate from a simple black 
body ( Hildebrand]| 1983] ). Although To is often set as a free pa¬ 
rameter, the models here lack sufficient photometry to suitably 
constrain this value. It was therefore assigned a fixed value of 
70pm, which is a reasonable value for such disks (Booth et al. 
|2013] >. A typical value of [3, also known as the spectral index, is 
two for ISM-like material. An equation of the form 


S v = LINko 




( 2 ) 


was therefore used, where Q is the solid angle of the emitting 
region, N the column density of dust, kq the opacity of the dust, 
and B.\(T) is the emission of a black body at a temperature T. 
With values of N, T, kq and spectral index (3 all being allowed to 
vary, the results become degenerate and therefore cannot be used 
to infer any specific disk properties. The total disk flux density 
however, is well approximated by this simple prescription. The 
stellar flux density was likewise fit assuming a pure black body 
spectrum of the star’s temperature and normalised to the J band 
flux. Both flux densities were then integrated according to Simp¬ 
son’s rule and divided to obtain the fractional disk luminosity. 

When disk excess is only detected at one wavelength, the 
value of Ld /L is extremely unconstrained. However, an upper 
limit was calculated by assuming a value of [3-2 and fitting the 
black body to the upper limits from ALMA and/or Herschel. 


3.5. SEDs and individual system parameters 
Systems with detections (2, 14, 17, 23) 

The SEDs for the four systems in which 1.3mm emission was 
detected are shown in Fig. |T] System 2 has a fractional disk lu¬ 
minosity that is very low, so is probably a massive debris disk. 
For the other three systems, however, the mid- to far-IR slope 
is suggestive of extended dust disks and not the thin belts often 
seen in debris disks. These three systems are therefore the most 
likely to be undergoing photoevaporation, and further observa¬ 
tions with a deeper detection threshold for the detection of gas 
will be invaluable. 

An interesting feature of these systems is that they are all 
believed to be in wide binaries, although there is only single- 
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epoch data available for these systems, so their binary nature 
has not been confirmed via proper motion (Kohler et al.|[2008 


Merfn et al.|2008 ; NguyeiTet al.|2012]|. Assuming that they are 


indeed binaries, one has to be careful in interpreting where the 
excess emission originates. In the case of System 2, the separa¬ 
tion is very large (>1400 au), and all detected emission can be 
confidently attributed to the one star. For System 17, the pro¬ 
jected separation is large enough that the photometry short of 10 
pm is resolved (the Spitzer resolution at 8 pm is ~2”), as is the 
mm photometry reported in this paper. The 24 and 70 pm points 
should be treated with some caution, but since no 1.3 mm emis¬ 
sion was detected for this alleged binary companion (see Fig. |2j, 
it is likely that emission at all wavelengths originates from a sin¬ 
gle circumstellar disk. A similar argument applies to System 14, 
although the ALMA image does suggest some kind of emission 
from the approximate location of the binary companion, which 
may contribute to the 24 pm flux found for this target as well. 
Given the wide separations of these targets, with projected sepa¬ 
rations of -420 and 560 au, it is unlikely that the tentative binary 
companion would influence the circumstellar disk evolution con¬ 
siderably. This is supported by Harris et al. ( |2012j >. who found 
no difference in disk luminosity for disks in wide binaries (sep¬ 
arations >300 au), compared to those around single systems. 

For System 23, the projected separation is too small for any 
of the detected excess emission to be resolved, and it could there¬ 
fore originate in either circumstellar disks or a single circumbi- 
nary disk. In either case, if the true separation of this system is 
<40 au, it is likely that any disks will undergo a very different 
evolution compared to disks around single stars. For example, 
observations suggest that binaries of separations <40 au inhibit 
the formation of protoplanetary disks (Kraus et al. 2012 1 . How¬ 
ever, disk-possessing close binaries can be found, and there is ev¬ 
idence to suggest that the tidal torque from such a binary could 
even slow down disk evolution ( | Alexander||20 \2) . This would 
cause prolonged lifetimes for the disk material, which may ex¬ 
plain this particular detection. However, the binary nature for all 
systems needs to be confirmed and the sample size increased be¬ 
fore conclusions can be drawn about how binarity affects this 
final stage of a disk’s lifetime. 

Systems 14 and 17 are the only systems in our sample that 
contain excess emission at 12 pm, suggesting the presence of 
dust located at small radii. Indeed, the other mid-IR detections 
of these two systems only deviate slightly from the average disk 
emission in Taurus. This, coupled with the fact that these two 
systems have the highest detected dust mass, makes them the 
most likely to have recently ceased accreting and be in the final 
stages of gas clearing. As mentioned in Section 3.2, the lack of 
a CO detection in these systems suggests that if this is the case, 
the gas-to-dust ratio must be severely depleted. 


Systems with large mid-IR excess (7, 9, 13, 21) 

Systems, 7, 9, 13, and 21 all have published 70 pm detection 
and steep rises in the mid-IR, which suggests a large population 
of cold dust (see Fig. [3]). As a result, the ALMA non-detections 
are rather surprising. We therefore include the black-body fits 
used to determine their fractional disk luminosity into Fig. 3] 
to investigate whether the value for the spectral index must be 
/3 »2. All fall below or on the ALMA upper limit with a spectral 
index [3- 2, with the exception of system 13, for which a value 
of [3 > 3.2 was required to satisfy the upper limits. This rather 
high value of beta suggests that the far-IR points of this system 
are not associated with the source, and |Matra et al.| ( [20121 ) show 
that show that emission at wavelengths longer than 70 pm could 
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Fig. 2. Images generated by the CLEAN algorithm for the four systems in which we detected continuum emission. The black cross denotes the 
position of the target, and the white cross the approximate position of any potential binary companions. The white ellipse denotes the beam size. 
The panel above each image displays spectra centred on the 12 CO(2-l) line for the region of continuum emission, which can be seen to contain 
only noise. 


be dominated by a source south-west of the system. Accounting 
for this, we only find a single black body of T=85K, and /3=2 is 
required to fit the emission of this system. 

For System 7, our age analysis found no PMS tracks that 
were consistent with its luminosity, suggesting that it is perhaps 
not a member of the young cloud and is already on its main- 
sequence. Massarotti et al. ( 2005j i support this hypothesis by 
showing that the proper motion of this system is indeed too high 
to be part of Taurus. The high proper motion makes it likely that 
this system is, in fact, a foreground main sequence star of a later 
spectral type with less extinction. 

Systems 9 and 21 both have excesses that begin at 70 gm, 
suggesting a population of dust far from the star. System 9 can 
easily be explained by material in a thin belt and is therefore 
probably a cold debris disk, albeit a very young one. The 160 n m 
point obtained for system 21 requires an even cooler debris disk 


of ~20K to explain all flux, corresponding to a distance of ~600 
au. This is extremely far from the star, and so this is unlikely to 
occur. Instead, it is more likely that the 160 gm point includes 
contamination from extended emission in the Herschel image 
from which it was derived (Cieza et al. 2012 1 . In any case, both 
systems fit with a spectral index of /3-2, which is consistent with 
dust similar to the ISM. 


If the long wavelength emission in Systems 13 and 21 is in¬ 
deed due to another source, then it is likely that all four of these 
systems are in their debris phase. This is also apparent from their 
fractional disk luminosities, which are approximately 10 5 . A 
deeper search with ALMA would clarify the situation for these 
uncertain objects since it could both detect and resolve the con¬ 
tinuum emission to a much smaller area than is possible with 
Herschel. System 7 is perhaps the least certain, however, with a 
relatively high fractional disk luminosity and some uncertainty 
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Log wavelength (//in) 


Fig. 3. SEDs for the four systems in our sample that display considerable mid-IR excess, but for which there was no 1.3mm detection. Symbols 
are identical to Fig. [I] with the addition of modified black body fits to determine if the emission is still consistent with one or two black bodies, 
taking the ALMA upper limits into consideration. 


on its spectral type. If it transpires that it is a late-type star, then 
the fractional disk luminosity will be increased, and this object 
would perhaps need to be re-classified. 


Likely debris disks 


For the remaining systems (Fig. [4]), only one or two excess points 
exist, and the upper limit on the 1.3mm flux is not so restric¬ 
tive. As a result it is easy to explain the excess as originating 
from a single temperature black body. The value of L^/L* is 
confined to the debris regime, but it is difficult to class these ob¬ 
jects as debris disks without knowledge of their gas content. The 
argument presented in section 3.2 suggests that, for the major¬ 
ity of disks, CO should have been detected despite the effects 
of photo-dissociation and our CO non-detections therefore sug¬ 
gest these disks are indeed gas poor. In support of this, previous 
studies have searched directly for FL tracers in similarly young, 
low mass systems and found nothing within a few au of the star 
(Pascucci et al. 2006} |Ingleby et al. 2009}. For example, |Pas- 


|cucci et ak ( |2006 1 were able to rule out gas masses above 0.04 
Mj within 3 au of the inner disk radius for the disks in their 
sample, which includes disks in a similar 5-15 Myr range. They 
conclude that gas dissipation is very efficient, so it would be sur¬ 
prising if the gas has survived in the low mass systems in our 
sample. We therefore consider it most likely that the majority 
of these systems are in their debris stage, containing only thin 
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belts of dust and no gas. Systems 1, 4, and 6, which all appeared 
too under-luminous to fit a PMS track, are probably also debris 
disks around background main-sequence stars. There is evidence 
to suggest that 70 pm observations are more sensitive to debris 
disks around such high-temperature objects (Cieza et al. 2008a), 
allowing them to be detected at larger distances. 


4. Discussion 

4.1. WTTS as young debris disks 

Figure [5] compares the ages and dust masses of the systems in 
our sample to those in the literature. Literature values for dust 
mass were all calculated from optically thin emission in a similar 
manner to this paper, and ages were estimated using either evo¬ 
lutionary track models (|Greaves et al.|2004| |Sheret et al.|2004| 
[Matthews et al.||2007[ |Ricci et ak 2010[ |Romero et al.||2012 1 , 


et al. 

2004 

Sheret et al. 2004;|Najita & Williams 2005;|Lestrade 

et al. 

2006 

), or cluster membership (Williams & Andrews 2006). 


The upper limits and detections for our sample clearly show that 
the dust masses for 23 of the 24 objects lie in the debris disk 
regime. Furthermore, the fractional disk luminosity for the ma¬ 
jority of the objects has a value of Lp/L* < 3 x 1CL 3 , suggesting 
they are most probably young debris disks (see Table [3}. Even 
System 2, in which there is a detection at 1.3 mm, has a value that 
places it in the debris regime, making it potentially a high-mass 
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Log wavelength (yu.m) 

Fig. 4. SEDs for the systems in our sample in which there was no 1.3mm detection, and for which the emission is consistent with a thin belt of 
material at constant temperature. Symbols are the same as in Fig. |T] 


debris disk rather than a photoevaporating disk. Systems 14, 17, 
and 23, on the other hand, have relatively high fractional disk 
luminosities, which, combined with their 1.3 mm detections, 
suggests they are not traditional debris disks and may instead 
be photoevaporating disks. The lack of gas detection in these 
systems may also be explained by photoevaporation, although 
another possibility worth bearing in mind is that they already 
possessed a reduced gas-to-dust ratio before reaching the WTTS 
state. Evidence for this possibility has recently been found by 
|Ansdell et al,| ( |2015) , in which a reduced gas-to-dust ratio of <2 
was found in an accreting system, potentially as a result of dust 
filtration ( |Rice et al.|2006| l. 


Taking the broad definition of a debris disk, where they are 
defined as gas-poor, geometrically thin dust disks at a uniform 
temperature, then the previously youngest known debris disks 
have ages in the range of 3-5 Myr ( Pascucci et al.|2006| |Chen| 
|et al.||2014| |Rigliaco et al.|[20l5] l. The systems observed here, 
however, generally have younger ages, with 10 systems lying 
within the age range 1.8 +qj to 4.7*0® Myr. Even accounting 
for the intrinsic uncertainty in these ages, they are significantly 
lower than previously detected debris disks, and this raises ques¬ 
tions as to why such young disks have not been observed be¬ 
fore. This is quite possibly the result of observational bias, since 
debris-disk emission is intrinsically faint and, as such, observa¬ 
tions are limited to nearby stars. The arrival of ALMA, however, 
now allows us to open this parameter space to young clusters at 
much greater distances of > 100 pc and our results clearly hint 
that rapid evolution into the debris phase is possible even for 


late spectral types. If confirmed, then these systems will lower 
the minimum age of debris disks and indicate a wider range of 
time scales for protoplanetary disk evolution. 

A more restrictive definition of debris disks describes them 
as containing second-generation dust, which is formed through 
a continuous process of collisions between planetesimals, and 
subsequent removal through radiation pressure and Poynting- 
Robertson drag ( [Wyatt[2008j >. Using this definition, it is unclear 
whether the disks in our sample can be classed as debris. Wah- 
haj et al.j( 2010j > suggest that dust around WTTS is primordial, 
based on the lower fraction of IR-WTTS found when comparing 
off-cloud to on-cloud sources. They suggest that this is caused 
by dissipating primordial dust since there is a link between sep¬ 
aration from the cloud and age, with the older WTTS located at 
increased separation from their parental cloud. This interpreta¬ 
tion is still rather uncertain however, because both on- and off- 
cloud sources in the survey have a wide variety of ages with 
considerable overlap, resulting in only a weak trend in age as 
a function of separation. Instead, the majority of both observa¬ 
tional evidence and photoevaporation models indicate that dust 
in the inner and outer regions of disks dissipates more or less si¬ 
multaneously (e.g. |Andrews & Williams|2005| |Alexander & Ar- 
mitage 2007), and most models of photoevaporation do not allow 
a considerable amount of dust to be left behind. This therefore 
implies that the debris disks seen around WTTS contain second 
generation dust. It has also been suggested that a large portion 
of the dust in protoplanetary disks could even be second gener¬ 
ation, based on the observation that the growth of grains from 
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Fig. 5. Comparison of the dust mass in our systems compared to those found in the literature. Cyan symbols denote the systems in our sample 
with detected 1.3 mm flux (numbered for clarity) and red symbols denote those those with upper-limits. Black symbols that lie above the dotted 


line correspond to known protoplanetary disks (data taken from 

Wyatt et al.|2003;; Ricci et al.|2010 

Romero et al. 2012). Black symbols in the 

lower portion correspond to known debris disks (data taken from 

Greaves et al. 2004, Liu et al. 2004 

Sheret et al. 2004;|Najita & Williams|2005 


|Lestrade et al.|20 06, Williams & Andrews 2006, Matthews et al. 2007). The horizontal dashed line is used to highlight the difference in dust mass 

between the populations 


micron to metre sizes can occur very rapidly (Dominik et al. 
2007[ >. This would quickly reduce the dust mass of protoplane¬ 
tary disks inferred from IR observations, and yet this value re¬ 
mains fairly constant for a range of disk ages ( jNatta et al.|2007| l. 
A mechanism of dust replenishment is therefore required in pro¬ 
toplanetary disks ( jDullemond & Dominik||2005j l, and we con¬ 
sider it likely that the debris disks in our sample contain second- 
generation dust as well. 

The observed fraction of the WTTS population that display 
an excess (~ 20%) is similar to the fraction of debris disks found 
around young FGK type stars (between 10-16%; [Hillenbrand 


|et al.| 2008 |Trilling et al.||2008| ). Although it is often assumed 
that IR-WTTS evolve into those without an excess, an alternative 
possibility is that these IR-WTTS mostly contain debris disks of 
second-generation origin, and these debris disks can then per¬ 
sist into their main sequence lifetime. The WTTS without ex¬ 
cess would then make up an entirely different population, and the 
difference between these two populations could be that the for¬ 
mer has a method of stirring the disk, which is a requirement for 
the collisional cascades that form debris disks. Stirring via stel¬ 
lar flybys (Kenyon & Broml ey|2002|) and self-stir ring via 1000 
km-sized planetesimals ( Kenyon & Bromley |2010 1 are both pos¬ 
sible causes for this. Perhaps the most exciting explanation for 
WTTS with debris disks, however, is that these systems have 
formed giant planets capable of stirring the disk (Mustill & Wy 


|att|2009| l. This would require giant planets orbiting at a few au, 
and the occurrence for such bodies around FGK stars has been 
estimated from radial velocity surveys at between 12% to 22% 
(|Lineweaver & Grether|2003]|Marcy et al.|2005]|Cumming et al.| 

2008). The similarity between this occurrence rate and the num¬ 
ber of WTTS systems that possess debris disks makes such sys¬ 
tems excellent places to perform planet searches. 
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4.2. Implications for photoevaporation models 

Some models of X-ray photoevaporation predict a sample of 
’relic’ transition disks with large cavities and high dust masses 
that persist for >10 Myr ( |Owen et al.[201 1| >. Our findings add to 
the growing body of research going against this prediction (e.g. 
|Cieza et al.|2010|[2012[|Mathews et al.|2012[ >, since the vast ma¬ 
jority of our IR-WTTS have no detectable 1.3 mm emission de¬ 
spite their ages being under 10 Myr. Some of the more recent X- 
ray photoevaporation theories introduce a mechanism they call 
"thermal sweeping", which disperses the remaining material in 
these massive dust disks in a small fraction (1-3%) of the disk’s 
total lifetime (Owen et al. |2013| > and removes the prediction of 
these relic disks. Likewise, the EUV models of photoevapora¬ 
tion predict short time scales for final disk clearing. EUV mod¬ 
els predict much lower photoevaporation rates in general, but 
are also only capable of forming transition disks when viscous 
accretion has already cleared a large amount of material. The 
resulting disks predicted by EUV models therefore have lower 
mass than in the X-ray case, and clearing can progress on time 
scales of between 1% to 10% of the disk’s lifetime ( Alexander! 
|& Armitage|2007| l. 

Previous surveys of WTTS have found that ~ 20% of WTTS 
display an excess that, when compared to the number of CTTS 
in the same regions, suggests that the disks around WTTS per¬ 
sist for 10% to 20% of the disk lifetime before moving into a 
diskless state ( |Cieza et al.|2007[|Wahhaj et al.|20fP) . This per¬ 
centage is somewhat higher than predicted by photoevaporation 
models, but its derivation assumes that all IR-WTTS follow the 
same evolution, moving from a CTTS, to an IR-WTTS, and fi¬ 
nally to a diskless state. If, as outlined above, the WTTS with 
debris disks are not in transition, then the apparent rarity of pho- 
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toevaporating disks will significantly lower this percentage and 
may bring it more in line with photoevaporation models. The 
small number of detections and lack of gas confirmation in our 
potential photoevaporating disks does not allow accurate esti¬ 
mation of this corrected percentage, but a wider survey of these 
photoevaporating disks may confirm this tentative result. 


4.3. Comparison to known debris disks with gas 


In recent years, a small population of debris disks have been 
found with detectable gas, leading to some debate as to its ori¬ 
gin. As with the dust, the origin of this gas is believed to be 
either primordial or formed through collisions of icy comet-like 
objects. One such system is 49 Ceti, with spectral type Al, a dust 
mass of ~0.3 M ffl and a 12 CO(2-l) integrated intensity of 2.0 Jy 
kms(Hughes et al. 2008). Its age has proved difficult to de¬ 
termine since it is not obviously a member of any associations, 
but |Thi et all] ( |200l[ ) believe it is a PMS star with an age of 8 
Myr, opening up the possibility that it is a high-mass analogy to 
the stars in our sample. If 49 Ceti were at the distance of Taurus, 
however, we would expect to detect this level of CO. Likewise, 
the 30 Myr-old, A4 type system HD 21997 is classed as a debris 
disk, containing only 0.09 M@ of dust, and yet it displays CO 
emission that we would have conclusively detected in our sur¬ 
vey ( jKospal et al.||2013| . If these systems are truly harbouring 
primordial gas, then the evolution for A-type stars must be dras¬ 
tically different to late-type stars to allow them to retain such a 
large quantity of gas at such low dust masses. Alternatively, the 
gas is secondary and both A type stars and most of the late-type 
stars in this survey lose their primordial gas by ages of <10Myr. 


5. Conclusion 

All the above sources in our sample are beyond the stage of ac¬ 
tive gas accretion. Their SEDs are suggestive of a depleted dust 
mass, and here we confirm this, either with their low 1.3 mm 
flux or with their non-detection from ALMA. The dust for all 
systems must, therefore, have either been largely removed or ag¬ 
glomerated into larger particles. The non-detection of CO lines 
in all systems is suggestive of a similar fate for the gas, which 
was probably removed by photoevaporation. Although photo¬ 
dissociation models for the disks studied here suggest that the 
CO abundance will be lowered slightly, for the majority of disks 
CO should have remained detectable. It is therefore probable that 
the gas-to-dust ratio has evolved in the majority of these disks to 
a value well below that of the ISM. The non-detection of CO 
is particularly surprising for Systems 14 and 17, whose SEDs 
strongly suggest large dust reservoirs. The depleted gas-to-dust 
ratio in these two systems may therefore be the result of a late 
stage of photoevaporation, or else they already possessed a re¬ 
duced gas-to-dust ratio before reaching the WTTS state. 

For those systems in which there was no 1.3 mm detection, 
it is probable that they are free from gas and contain dust masses 
and distributions similar to debris disks. This is apparent from a 
comparison of their fractional disk luminosities and dust mass 
upper-limits to that of known debris disks, as both lie in the 
debris disk regime. These systems, however, are clearly much 
younger than the majority of debris disks allowing for more strict 
constraints on debris disk formation time scales than ever before. 
A deeper study with ALMA will be invaluable to determine con¬ 
clusively their evolutionary state, as well as to confirm the dust 
masses of the suspected debris disks. 
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